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Abstract We use our Monte Carlo radiative transfer code to study non-embedded prestel- 
lar cores and cores that are embedded at the centre of a molecular cloud. Our 
study indicates that the temperature inside embedded cores is lower than in iso- 
lated non-embedded cores, and generally less than 12 K, even when the cores are 
surrounded by an ambient cloud of small visual extinction (Ay ~ 5). Our study 
shows that the best wavelength region to observe embedded cores is between 
400 and 500 jim, where the core is quite distinct from the background. We also 
predict that very sensitive observations (~ 1 — 3 MJy sr~') at 170-200 yum can be 
used to estimate how deeply a core is embedded in its parent molecular cloud. 
Finally, we present preliminary results of asymmetric models of non-embedded 
cores. 



Introduction 

Prestellar cores are cores that are either on the verge of collapse or already 
collapsing (e.g. Myers & Benson 1983, Ward-Thompson et al. 2002). They 
represent the initial stage of star formation and their study is important since 
theoretical models of star formation are very sensitive to the initial conditions. 

Prestellar cores have been observed either isolated or embedded in pro- 
toclusters. Isolated prestellar cores (e.g. L1544, L43, L63) have extent > 
1.5 x 10 4 AU and masses 0.5 - 35 M© (Ward-Thompson et al. 1999, Andre et 
al. 2000). On the other hand, prestellar cores embedded in protoclusters (e.g. 
in p Oph, NGC2068/2071) are generally smaller, with extent ~ 2 -4 x 10 3 AU 
and masses ~ 0.05 - 3 M© (Motte et al. 1998, Motte et al. 2001). 
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In this paper, we present radiative transfer models of non-embedded and em- 
bedded prestellar cores, performed using a 3-D Monte Carlo radiative transfer 
code we have developed (Phaethon). 

1. The Method: Monte Carlo Radiative Transfer 

Our method (Stamatellos & Whitworth 2003) is similar to that developed 
by Wolf, Henning & Stecklum (1999) and Bjorkman & Wood (2001). We 
represent the radiation field of a source (star or background radiation) by a large 
number of monochromatic luminosity packets (L-packets). These L-packets 
are injected into the system and interact stochastically with it. If an L-packet is 
absorbed its energy is added to the local cell and raises the local temperature. 
To ensure radiative equilibrium the L-packet is re-emitted immediately with 
a new frequency chosen from the difference between the local cell emissivity 
before and after the absorption of the packet (Bjorkman & Wood 2001). This 
method conserves energy exactly, accounts for the diffuse radiation field and its 
3-dimensional nature makes it attractive for application in a variety of systems. 

The code has been thoroughly tested using the thermodynamic equilibrium 
test (Stamatellos & Whitworth 2003) and also against benchmark (Ivezic et al. 
1997) and previous (Bjorkman & Wood 2001) calculations. 

2. Non-Embedded Prestellar Cores 

We represent prestellar cores by Bonnor-Ebert (BE) spheres, in which grav- 
ity is balanced by gas pressure (Bonnor 1956, Ebert 1955). In many cases 
this is a good approximation to prestellar cores (e.g. Alves et al. 2001, Ward- 
Thompson et al. 2002). 

For the radiation incident on the core, we use the Black (1994) interstellar 
radiation field (hereafter BISRF), that consists of an optical component due to 
radiation from giant stars and dwarfs, a component due to thermal emission 
from dust grains, mid-infrared radiation from non-thermally heated grains and 
the cosmic background radiation. 

The opacity of the dust at the low temperatures (5-20K) and high densities 
(10 4 - 10 7 cm~ 3 ) expected in prestellar cores, is quite uncertain. In our study, 
we use the opacities calculated by Ossenkopf and Henning (1994), for grains 
that have coagulated and accreted thin ice mantles. 

We find that the temperature inside non-embedded cores drops from around 
17 K, at the edge of the core, to a minimum at the centre, which maybe as 
low as 7 K, depending on the visual extinction to the centre of the core. Our 
results are similar to those of previous studies (Evans et al. 2001, Zucconi et 
al. 2001). 
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3. Embedded Prestellar Cores 

The radiation field incident on cores embedded in molecular clouds is differ- 
ent from that incident on isolated non-embedded cores. The ambient molecular 
cloud absorbs the UV, optical and NIR part of the radiation and reemits it in 
the FIR and submm region (Mathis et al. 1983). It also makes the radiation 
incident on the core anisotropic, because, in general, the ambient cloud is not 
homogeneous. Another factor that contributes to the anisotropy of the radia- 
tion incident on an embedded core is the presence of stars or protostars in the 
vicinity of the core or the cloud (e.g as in p Oph; Liseau et al. 1999). 

The Model. In this first approach (Stamatellos & Whitworth 2003), we 
study a spherical core at the centre of a molecular cloud of uniform density. 
The radiation incident on the molecular cloud is the BISRF, but the radiation 
incident on the core is enhanced in the FIR and submm and reduced at shorter 
wavelengths, as a result of the presence of the molecular cloud around the core. 

We chose the parameters of our models so as to mimic the embedded prestel- 
lar cores and the conditions in the p Oph protocluster (see Motte et al. 1998): 
core sizes 4 - 8 x 10 3 AU, masses 0.4 - 1.2 M©, ambient cloud particle density 
« tot = 0.96 x 10 4 cm~ 3 and ambient cloud pressure ~ 10 6 cnT 3 K. In Figs. 1-2, 
we present our calculations for a supercritical core at gas temperature T-15 K 
with mass 0.8 M©, under external pressure P ext - 10 6 cm -3 K, surrounded by 
a spherical ambient cloud with different visual optical depths. 

Core Temperature Profiles. We find that the presence of even a mod- 
erately thick cloud (Ay = 5) around the core, results in a less steep temperature 
profile inside the core than in the case of a core that is directly exposed to the 
BISRF. When there is no surrounding cloud (Fig. 1, dashed lines), the temper- 
ature drops from -16 K at the edge of the core to around 6-7 K in the centre 
(AT « 9 - 10 K, depending on the core density), whereas with a Ty = 5 ambi- 
ent cloud (Fig. 1 , dotted lines) the temperature drops from around 11 K to 7 K 
(AT « 4 K). Our studies show that dust temperatures inside embedded cores 
are probably lower than 12 K in cores surrounded by even a relatively thin 
cloud (Ay ~ 5), which seems to be the case for many of the prestellar cores 
in p Ophiuchi. Previous studies (Motte et al. 1998, Johnstone et al. 2000) of 
cores in this region, assumed isothermal dust at temperatures from 12 to 20 K, 
when calculating core masses from mm observations and, thus, they may have 
underestimated masses by a factor of 2. 

Recent studies (Andre et al. 2003, Andre, this volume) find similar temper- 
ature profiles, using a different approach, in which they estimate the effective 
radiation field incident on an embedded core from observations. 
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Figure 1. (a) Density, (b) dust temperature, and (c) SEDs, for a supercritical core (see text) 
surrounded by a spherical cloud with visual optical depth 20 (solid lines), 5 (dotted lines) and 
(dashed lines). The dash-dot line on the SED graph corresponds to the background SED. 



SEDs and Intensity Profiles. At 90 microns the core is seen in ab- 
sorption against the background (Fig. 2a), and the intensity increases towards 
the edge of the core. For very centrally-condensed cores the decrease towards 
the centre is just -8-10 MJy sr _1 , but for less centrally-condensed cores it 
is even lower. Thus, very sensitive (say - 1 - 3 MJysr -1 ) observations are 
needed to detect cores in absorption at 90 jum. 




Figure 2. Intensity profiles at (a) 90, (b) 170, and (c) 450 and 850 /im, for the models in 
Fig. 1. The horizontal solid lines on the profiles correspond to the background intensity at the 
wavelength marked on the graph. 

At wavelengths near the peak of the core emission (150-250 /mi) the inten- 
sity increases by a small amount (-5-20 MJysr -1 above the background) 
towards the edge of the core (Fig. 2b). The higher the increase in the inten- 
sity near the core boundary, the less embedded is the core. Thus, very sensitive 
observations of embedded prestellar cores at 170-200 /mi, might allow us to de- 
termine the extinction of the cloud surrounding the core, and, thus, to estimate 
roughly the position of the core inside the cloud. However, more sophisticated 
modelling is required, with an accurate density profile for the cloud and taking 
into account the close environment of the core under study. 
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Finally, at submillimeter and millimetre wavelengths (400-1300 /im) the 
intensity drops toward the edge of the core considerably (Fig. 2c). The core 
can be easily observed at 400-500 /mi, where the contrast with the background 
is quite large (~ 50 - 150 MJy sr _1 ). At wavelengths longer than ~ 600 yum the 
background radiation becomes important and the core emission is not much 
larger than the background emission (~ 20-50 MJy sr -1 larger). 

Diagnostics. In Table 1, we list the peak intensities at various wave- 
lengths for cores embedded in molecular clouds with visual optical depths 
5 and 20. This table indicates that embedded cores are most easily distin- 
guished from the background radiation around 450 yum. The peak emission 

Table 1. Typical peak* intensities for embedded cores 



A (pm) 


//(MJysr 1 ) 




f~cloud — 5 f cloud — 20 


90 b 


5-15 ~3 


170 


10-15 ~ 3 


450 


55-160 40-130 


850 


20-80 15-70 


1300 


10-40 10-25 



* The term peak refers to the maximum intensity above or below the background (as noted) at a specific 
wavelength. 

a Approximate peak intensities for a core embedded in a cloud with visual optical depth 5 and 20. The 
deeper the core is embedded the less distinct from the background is. The lower value corresponds to a 
subcritical core and the higher value to a supercritical (i.e. more centrally condensed) core. 
b At 90 p.m the core seen in absorption against the background. 

from embedded cores could be as low as ~ 10 MJy sr 1 above the background 
at 1300 yum, but it's at least ~ 40 MJysr -1 at 450 yum. The wavelength range 
between 400-500 seems favourable for observing embedded cores but the at- 
mospheric transmission is not good in this range and space observations are 
needed. The upcoming Herschel (to be launched in 2007) will be operating in 
this range. 

4. Asymmetric Models of Non-Embedded Cores 

The model. Generally, cores are not spherically symmetric. Here, we 
assume non-embedded cores with a disk-like asymmetry, i.e. the core is denser 
on the equatorial plane; 
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where p c is the the density at the centre of the core, ro the scale length, and 
A is a factor that determines how asymmetric is the core. This density profile 
(Fig. 3, left) resembles the BE sphere density profile (drops as r' 2 at large radii 
and gets flatter near the centre). In Fig. 4, we present our calculations for a core 
with n c = 10 6 cnT 3 , r = 2 x 10 3 AU, R core = 2 x 10 4 AU, M = 7.3 M Q and 
e = 2.5 (e is the ratio of the optical depth at 9 - 90, i.e. looking at an edge-on 
core, to the optical depth at 9 = 0). 




Figure 3. (a) Density, (b) dust temperature, and (c) SEDs, for an asymmetric core (see text). 
The core is divided into 15 equal-angle cells (each cell's angle span is 12°). The lower density 
cell corresponding to 6 = 0° (lower curve on the density plot), has higher temperature (upper 
curve in the temperature plot). The dotted line on the SED graph corresponds to the incident 
SED. 

Temperature Profiles. The dust temperature is 9 dependent (Fig. 3, 
centre). As expected, the 'equator' of the core is colder than the 'poles'. The 
difference in temperature is around 3-4 K for the core under study (e = 2.5). 
The difference is expected to be larger for less symmetric cores. 

SEDs and Images. The SED distribution (Fig. 3, right) of the core is the 
same at any viewing angle, because the core is optically thin to the radiation 
it emits (FIR and longer wavelengths). However, the isophotal maps are quite 
different for different viewing angles (Fig. 4). 

At 200 /j,m the outer, hotter parts of the core dominate the core emission. 
The core appears spherical when viewed pole-on and elongated when viewed 
edge-on. The image at 30° is quite interesting, with two absorption 'blobs' 
near the centre of the core. A quick search through the Kirk (2003) sample of 
ISO/ISOPHOT observations did not reveal cores with such distinctive features, 
which is surprising, because one excepts some of the observed cores to be 
viewed at such angles. Further study is required to see whether this is due to 
the sensitivity and resolution of observations, due to selection effects, or it is 
connected with the core structure and its environment (i.e. the ambient cloud). 
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Figure 4- Isophotal maps at 200 (left panel) and 850 /im (right panel) at viewing angles 
0°, 30° and 90° (top to bottom), for a disk-like asymmetric core (see text). The core appears 
elongated when viewed at a direction different from 6 = 0°. The colour image is available on 
the CD-ROM. 
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At longer wavelengths, like 850 /urn, the core emission is regulated by the 
column density. Thus, when the core is viewed edge-on the intensity is larger 
at the centre. For the same reason the core looks elongated when the observer 
is looking at it at any other direction than pole-on. 

5. Conclusions 

We applied a Monte Carlo radiative transfer method to spherical cores em- 
bedded inside molecular clouds. We find that the temperature inside these 
cores is less than 12 K, even for an ambient cloud with moderate visual extinc- 
tion ~ 5. We also studied asymmetric non-embedded cores and the preliminary 
results show that a small disk-like asymmetry in the density distribution will 
make the core look elongated when viewed at a random angle. 
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